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ABSTRACT 

ridally  induced,  sudden  temperature  drops  in  Nan  Wan,  partieularly  pronouneed  during  the  spring  tide,  often  bear 
eeologieal  eonsequenees.  We  u.se  an  Oeean  Noweast/Forecast  System  (ONFS)  at  the  US  Naval  Research  Laboratory  (NRL)  to 
examine  a  widely  publicized  cold-water  intrusion  event  in  July  2008,  and  compare  it  with  what  we  perceived  as  a  normal 
intrusion  event  in  June  2007.  For  a  normal  cold-water  intrusion,  cold  anomalies  eventually  propagate  away  in  forms  similar  to 
internal  Kelvin  waves.  For  early  July  of  2008,  the  NRL  ONFS  suggested  an  anomalous,  subtidal  and  es.sentially  eastward 
current  near  Nan  Wan.  The  subtidal  current,  superimposed  on  the  ebb  current,  blocks  the  propagation  of  cold  anomalies  in  the 
direction  of  internal  Kelvin  waves,  inducing  anomalous  and  ineftlcient  cold-water  dispersal.  The  NRL  ONFS  also  suggested 
the  anomalous  subtidal  current  as  a  post -typhoon  response  to  the  passage  of  Typhoon  Fengshen. 

Key  word.s:  Nan  Wan  upwelling.  Tide- induced  upwelling.  Cold-w  ater  intni.sion,  Typhoon-Ocean  interaction 

Ciiathm  Ko.  D.  S..  S.  V.  Chao.  P.  Huang,  and  S.  F.  Lin.  2009:  .Anomaloiis  upwelling  in  Nan  Wan:  Julv  200S.  Ten'.  Almos.  Ocean  Si  i..  20,  S39-S52.  doi: 
I0  3319/TA0  200S.  I  l.25.0l(Oc) 


1.  INTRODUCTION 

Nan  Wan,  or  the  south  bay  of  Taiwan,  is  a  semi-enclosed 
basin  open  to  the  Luzon  Strait  (Fig.  1 ).  The  northward  flow¬ 
ing  Kuroshio  follows  the  w  estern  extreme  of  the  Philippine 
Sea  from  the  east  coast  of  Luzon  to  the  east  coast  of  Taiwan. 
In  summer,  the  Luzon  Strait  is  more  or  le.ss  free  from  Ku- 
ro.shio  encroachment  and,  in  consequence,  the  Kuroshio 
often  stays  some  distance  away  from  Nan  Wan  and  its  vi¬ 
cinity.  In  winter,  the  encroaching  Kuroshio  often  enters  the 
northern  South  China  Sea,  looping  anticyclonically  before 
exiting  the  South  China  Sea  along  the  southern  tip  of  Taiwan 
(e.g.,  Centurioni  et  al.  2004;  Wu  and  Chiang  2007;  Liang  et 
al.  2008).  Conceivably,  the  encroachment  breeds  eddy  shed¬ 
ding  and  moves'  the  Kuroshio  much  closer  to  Nan  Wan  and 
its  vicinity. 

Tidal  currents  to  the  south  of  Nan  Wan  are  essentially 
eastward  during  ebb  and  westward  during  flood.  Subtidal 
currents,  especially  tho.se  related  to  the  Kuroshio,  could  con¬ 
ceivably  modify  the  .strength  of  the  flood  and  ebb  currents. 
Diurnal  tides  are,  more  often  than  not,  stronger  than  semi¬ 
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diurnal  tides  in  the  Luzon  Strait  and  Nan  Wan.  The  ano¬ 
malous  amplification  of  diurnal  tides  over  semidiurnal 
tides  is  due  to  the  basin-wide  resonance  that  involves  the 
entire  basin  of  the  South  China  Sea  (Jan  ct  al.  2007).  In  Nan 
Wan  and  its  vicinity,  the  diurnal  tide,  modulated  by  a  se¬ 
midiurnal  beat,  often  induces  upwelling  or  cold-water  in¬ 
trusion  that  occurs  daily  (Lee  et  al.  1997).  The  upwelling  is 
particularly  strong  and  widespread  during  the  spring  tide 
(Chen  et  al.  2004).  On  the  ea.st  side  of  Nan  Wan,  the  tem¬ 
perature  drop  lasts  for  a  few  hours  each  and  usually  occurs 
once  every  diurnal  tidal  cycle.  In  west  and  central  parts  of 
the  basin,  similar  temperature  drop  usually  occurs  tw  ice  in 
one  diurnal  tidal  cycle  (Lee  et  al.  1999).  The  temperature 
drop  could  be  sudden  and  drastic.  On  24  November  1988, 
for  example,  a  sudden  temperature  drop  of  more  than  1 0°C 
resulted  in  a  mass  kill  of  fish  (Su  et  al.  1 989). 

Theoretically,  both  capes  (Mou-Bi-Tou  and  O-Luan-Bi 
in  Fig.  1  b)  and  their  southward  extension  ridges  are  capable 
of  producing  upw  elling.  A  ridge  induces  upwelling  by  lift¬ 
ing  the  incoming  tidal  current.  A  cape  can  achieve  a  similar 
effect.  It  blocks  the  incoming  tidal  current,  which  needs  to 
detour  around  the  cape  to  supply  w^aters  to  the  lee  side.  Tlie 
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inefficiency  of  water  supply  creates  a  relatively  low  sea  level 
anomaly  in  the  lee  and  bottom  waters  normally  upwell  to  fill 
the  void.  This  is,  in  essence,  the  conclusion  from  the  pre¬ 
vious  numerical  model  (Lee  et  al.  1999).  Given  the  com¬ 
plex  tidal  constituents,  bathymetry  and  coastline  orientation 
in  and  around  Nan  Wan,  it  is  diillcult  to  assess  the  impor¬ 
tance  of  each  cape  and  ridge  in  producing  upwelling  without 
realistic  numerical  simulations. 

In  the  afternoon  of  3  July  2008  local  time  (early  morning 
of  July  3  in  GMT),  fishemien  spotted  a  large  school  of  fish 
floating  in  the  western  reaches  of  Nan  Wan.  The  collective 
diagnosis  afterwards  pointed  to  hypothermia  as  the  likely 
cause.  A  bottom  temperature  sensor  at  station  T  lent  further 
support  to  this  idea.  The  event  motivated  us  to  investigate 
further  using  an  Ocean  Nowcast/Forecast  System  (ONFS)  of 
the  US  Naval  Research  Laboratory  (NRL)  (Ko  et  al.  2008). 
Not  surprisingly,  the  NRL  ONFS  identified  this  period  as  a 
spring  tide  that  generated  strong  upwelling  in  Nan  Wan. 
When  comparing  the  July  2008  upwelling  event  with  pre¬ 
vious  ones,  wc  also  identified  two  types  of  pathways  for  the 
dispersal  of  cold  anomalies:  one  disperses  in  the  propagation 
direction  of  internal  Kelvin  waves  and  the  other  in  the  oppo¬ 
site  direction.  In  the  formal  case  or  nonnal  upwelling  event, 
cold  anomalies  .seem  to  disperse  easily  and  do  not  over¬ 
accumulate  in  Nan  Wan.  In  the  latter  ca.se  or  anomalous  up¬ 
welling  event,  the  dispersion  seems  slow  and  the  temp>era- 
ture  drop  in  Nan  Wan  becomes  more  severe.  For  the  par¬ 
ticular  cold  event  in  July  2008,  the  NRL  ONFS  suggests 
that  typhoon  Fengshen  in  the  preceding  month  created  the 
preconditioning  for  the  anomalous  upwelling. 

2.  OCEAN  NOWCAST/FORECAST  SYSTEM 

We  use  a  real-time  Ocean  Nowcast/Foreca.st  Sy.stem 
developed  at  the  US  Naval  Research  Laboratory  to  facili¬ 
tate  the  present  study.  The  NRL  ONFS  produces  a  daily, 
short-term  (under  a  week)  forecast  of:  mesoscale  ocean  cur¬ 
rents,  temperature,  salinity,  and  sea  level  variation.  The  sys¬ 
tem  can  be  relocated  to  ditTerenl  regions  and,  once  set  up  for 
a  particular  region,  operates  automatically.  The  system  is  an 
integration  of:  a  data-assimilating,  dynamical  ocean  model, 
a  statistical  data-analysis  model,  and  various  data  streams 
for  ocean  bathymetry,  climatological  data,  surface  forcing, 
open  boundary  forcing,  and  observations  for  data  assimila¬ 
tion  (Ko  et  al.  2008).  The  NRL  Modular  Ocean  Data  Assimi¬ 
lation  System  (MODAS;  Carnes  et  al.  1 996;  Fox  et  al.  2002) 
is  used  within  ONFS  as  the  data-analysis  model.  MODAS 
uses  satellite  data,  in-situ  obser\  ations  and  histoncal  statis¬ 
tics  to  generate  three-dimensional  ocean  temperature  and 
salinity  analyses.  The  analyses  are  then  assimilated  into  the 
dynamic  model  to  produce  an  ocean  nowcast.  From  the  now- 
cast,  the  forecast  is  conducted  without  data  assimilation  us¬ 
ing  a  meteorological  forecast.  This  paper  only  uses  the  now- 
cast  portion  of  the  results  for  analysis. 


This  paper  primarily  uses  results  from  a  high  resolution 
Luzon  Strait  NowcastTorecast  System  (LZSNFS),  an  appli¬ 
cation  of  NRL  ONFS,  to  discuss  upwelling  in  and  around 
Nan  Wan.  The  system  domain  extends  from  1 1 6®E  to  1 25®E 
and  from  17°N  to  24°N,  covering  the  northern  South  China 
Sea,  the  Luzon  Strait,  and  a  portion  of  the  Philippine  Sea 
with  a  horizontal  resolution  of  2.3  km.  In  terms  of  vertical 
resolution,  there  are  1 1  a-layers  in  the  top  147  m  and  29 
z- layers  below.  The  system  derives  its  ocean  bottom  topo¬ 
graphy  from  the  Digital  Bathymetry  Data  Base  with  2- 
minute  resolution  (DBDB2).  Open  boundary  conditions  are 
derived  from  a  larger  scale  ONFS  for  the  entire  East  Asian 
Seas  (EASNFS).  The  Oregon  State  University  (OSU)  tidal 
forcing  drives  the  tidal  currents.  A  global  weather  forecast 
model  (Navy  Operational  Global  Atmospheric  Prediction 
System  or  NOGAPS)  and  a  regional  weather  forecast  mo¬ 
del  (Coupled  Ocean/Atmosphere  Mesoscale  Prediction  Sys¬ 
tem  or  COAMPS)  provide  the  ocean  surtace  forcing.  The 
LZSNFS  also  assimilates  satellite  altimeter  data  and  Multi- 
Channel  Sea  Surface  Temperature  (MCSST)  from  satellite- 
derived  AVHRR  to  improve  accuracy  further. 

The  system  is  able  to  produce  rather  accurate  tidal  cur¬ 
rent  in  the  Luzon  Strait  and  northern  South  China  Sea.  For 
example,  tidal  currents  flowing  over  the  ridges  in  the  Luzon 
Strait  generate  large-amplitude  nonlinear  internal  waves  in 
the  northern  South  China  Sea,  vividly  seen  in  satellite  im¬ 
ages.  To  reproduce  the  timing  and  location  of  these  internal 
waves  requires  accurate  spatial  and  temporal  distribution  of 
tidal  currents  and  LZSNFS  was  able  to  achieve  this  (Chao  et 
al.  2007). 

In  addition  to  LZSNFS,  we  also  use  results  from  the 
larger  scale  EASNFS,  another  application  of  NRL  ONFS,  to 
diagnose  the  origin  of  anomalous  Nan  Wan  upwelling  in 
July  2008.  The  EASNFS  covers  the  western  Pacific  and  all 
the  Asian  marginal  seas  from  17.3°S  to  52.2®N  and  from 
99.2°E  to  1 58.2®E  with  a  resolution  of  1/1 6  degree  or  9.2  km 
at  latitude  of  Luzon  Strait.  The  system  derives  its  open  ocean 
boundary  conditions  from  NRL  global  NCOM  (Navy  Co¬ 
astal  Ocean  Model)  (Rhodes  et  al.  2002)  that  operates  daily. 
Contrary  to  fine-scale  LZSNFS,  the  EASNFS  does  not  si¬ 
mulate  tides.  The  vertical  resolution,  forcing  and  methodo¬ 
logy  of  EASNFS  are  otherwise  similar  to  that  in  LZSNFS. 
The  EASNFS  has  been  applied  for  studies  of  current  at 
Korea/Tsushima  Strait  (Teague  et  al.  2006),  transport  in  the 
Gulf  of  Papua  (Keen  et  al.  2006),  typhoon-ocean  interaction 
al  the  West  Pacific  Ocean  (Lin  et  al.  2008)  and  air-sea  inter¬ 
action  at  the  Japan/East  Sea  (Mooers  et  al.  2009). 

3.  NORMAL  UPWELLING  DURING  SPRING 
TIDES 

Figure  2  illustrates  the  agreement  between  the  observed 
sea  level  and  LZSNFS-produced  counterpart  at  station  T 
from  15  May  to  15  June  2007.  Being  a  mixture  of  diurnal 


Anomalous  Upwelling  in  Nan  Wan:  July  2008 


841 


laj  t2i  122  123  124  ^  ^^20  3S'  120  40  120  4&^  *20  50  *20  5S  12100 

22  06 


aaoo 


21  56 


2V50 


2r*5 

lia  110'  120  121  ■  122  123  12<  120"a5'  120*40-  »20*4S'  %20'5ff  IfOW  12100 


5000  4000  3000  2000  1000  0  1200  1000  000  4W0  -400  JOT  0 

B«ttiyrn«try  (m)  Bitiiv'netrv  (m} 

Fig.  ! ,  Bathymetry  in  the  vicinity  or(ii)  Luzon  Strait  and  (b)  Nan  Wan.  Two  capes  bounding  Nan  Wan  are  Mou-Bi-Tou  (MBT)  and  O-Luan-Bi  (OLB). 
Station  T  in  (h)  records  bottom  water  temperature  at  7  m  depth  and  sea  level  variations. 


t,2 
LO 

o.a 

0.« 

S  0-^ 

I 

I  0.0 

u 

>  -0.2 

I  '0.4 

g  -0^6 

1/1  -o.i 
-to 
-1.2 

05/15  06/01  06/15 

2007 

Fig.  2.  The  observed  sea  level  s  anation  (in  red)  and  model-produced  counterpart  (in  blue)  at  station  T  from  15  May  to  15  June  2007 


NanWan 


I  I  1  I  I  1  l  i~Y~i  I  I  I  I  I  I  I  I  I  I  I  i  I  I  1  i  i  I  i  I 


and  semidiurnal  tides,  the  sea  level  variation  also  eontains  a 
pronouneed  spring-neap  eyele.  We  ehoose  June  2  to  illus¬ 
trate  Nan  Wan  upwelling  event  during  the  spring  tide.  Tidal 
range  in  this  day  exeeeds  1 .6  m. 

Figure  3  shows  bihourly  eurrents  and  temperature  at 
25  in  depth  for  2  June  2007.  At  hour  0  (GMT),  the  .start  of 
ebb  eurrent,  upwelling  begins  to  appear  southwest  of  O- 
Liiaii-Bi.  The  southward  extension  ridge  from  0-Luan-Bi 
apparently  lifts  the  ebb  eurrent  to  produee  the  upwelling.  It 


intensifies  and  expands  until  hour  8.  Upwelling  off  the  tip  of 
Mou-Bi-Tou  .starts  at  hour  4,  and  thereafter  intensifies  and 
expands  until  hour  8.  At  hour  10,  as  the  tidal  euiTent  re- 
eesses,  the  eold  anomaly  off  Mou-Bi-Tou  begins  to  move  to 
the  western  eoa.st  of  Taiwan,  while  the  eorresponding  ano¬ 
maly  otTO-Luan-Bi  moves  to  the  eentral  basin  of  Nan  Wan. 
Both  movements  are  in  keeping  with  the  propagation  direc¬ 
tion  of  coa.stal-trapped  waves.  From  hour  10  to  hour  14,  the 
two  cold  anomalies  continue  to  propagate  away  in  the  pro- 
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pagation  direction  of  coastal-trapped  waves.  By  hour  16, 
cold-water  masses  have  completely  moved  out  of  Nan  Wan. 
fhc  right-bounded  movement  of  the  two  cold-water  masses 
explains  why  the  temperature  drop  occurs  once  daily  on  the 
east  side  of  Nan  Wan.  but  twice  daily  in  the  west  and  central 
parts  of  the  basin  (Lee  et  al.  1999).  Leaving  Nan  Wan  up- 
welling  aside,  the  ebb  current  also  generates  a  cold  cyclonic 
eddy  on  the  lee  (cast)  side  of  0-Luan-Bi.  It  begins  to  form  at 
hour  4,  expands  and  intensifies  until  hour  12,  and  decays 
slowly  therealter. 

In  coastal  oceans  with  continuous  stratification  and 
variable  bottom  slope,  eoastal-trapped  waves  become  a  hy¬ 
brid  of  internal  Kelvin  weaves,  continental  shelf  waves  and 
topographic  Rossby  waves  (Wang  and  Mooers  1976).  It  is 
dilTieult  to  separate  them  cleanly  even  if  the  oceanographic 
.setting  is  strictly  unifomi  alongshore.  Ilie  irregular  coastline 
and  bathymetry  around  Nan  Wan  complicate  the  setting 


ev  en  further.  Nevertheless,  further  results  point  to  the  close 
resemblance  of  cold-anomaly  propagation  to  internal  Kelvin 
waves.  In  the  LZSNFS,  sizable  cold  anomalies  and  their 
right-bounded  propagation  show'  up  together  during  the 
spring  tide  but  diminish  thereafter.  Tlie  predominance  of 
density  forcing  associates  the  right-bounded  propagation 
mostly  to  internal  Kelvin  waves,  which  owe  their  existence 
to  density  anomalies,  rather  than  continental  shelf  waves  or 
topographic  Rossby  waves,  which  arise  from  the  bottom 
slope.  Judging  empirically,  one  can  estimate  the  phase  speed 
from  the  movement  of  the  cold  anomaly  that  initially  ap¬ 
pears  off  Mou-Bi-Tou.  From  hour  8  to  hour  14,  the  propaga¬ 
tion  speed  is  about  45  km  per  day  or  --0.5  m  s'*.  Ambient 
currents  around  Nan  Wan,  tidal  and  .subtidal  alike  do  not 
seem  able  to  stop  the  right-bounded  propagation.  To  the  con¬ 
trary,  the  right-bounded  propagation  thrives  in  the  “shadow 
zone"  of  the  adverse  ambient  currents.  For  example,  from 


Fig  Bihourly  currents  and  temperature  at  25  m  depth  for  2  June  20()7,  showing  nonnal  cold-water  intrusion  during  the  spring  tide.  Numbers  on  land 
indicate  time  (GMT)  in  hours. 
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liour  6  to  hour  10,  the  westward  expansion  of  the  eold  ano¬ 
maly  from  O-Luan-Bi  is  restricted  mostly  to  Nan  Wan  and 
its  immediate  south,  where  the  adverse  ambient  current  is 
weak  or  reversed.  From  hour  12  to  hour  16,  the  right- 
bounded  propagation  from  Mao-Bi-Tou  to  the  west  coast  of 
Taiwan  occurs  at  the  expense  of  diminishing  adverse  am¬ 
bient  eurrent. 

The  right-bounded  propagation  of  cold  anomalies  in 
Fig.  3  is  against  the  prevailing  tidal  eurrent  from  hour  6  to 
hour  10.  Tidal  current  is  therefore  not  responsible  for  it.  To 
associate  it  with  the  earth’s  rotational  effeet,  we  have  termi¬ 
nated  the  Coriolis  force  in  Fig.  4  after  hour  6,  the  onset  of 
right-bounded  propagation.  By  removing  the  earth  rotation, 
tidal  currents  also  change  gradually  in  time.  Nevertheless, 
the  right-bounded  propagation  diminishes  quickly  before 
the  tidal  eurrent  changes  signifieantly.  For  example,  the  tidal 
eurrent  in  Fig.  4  (without  the  earth  rotation)  essentially  re¬ 


mains  the  same  as  in  Fig.  3  (with  the  earth  rotation)  from 
hour  6  to  hour  1 2,  but  the  right-bounded  propagation  essen¬ 
tially  vanishes  in  Fig.  4  within  only  six  hours  after  the  termi¬ 
nation  of  Coriolis  force.  This  is  plausible,  beeause  the  earth 
rotation  supports  the  propagation  of  internal  Kelvin  waves. 

If  idealized  to  the  extreme,  internal  Kelvin  waves  pro¬ 
pagate  along  a  vertical  coast  and  over  a  flat  bottom.  With  a 
seaward  deepening  bottom  in  the  present  .setting,  where  is 
the  effeet ive  eoastal  wall  and  appropriate  water  depths  to 
support  the  cold-water  propagation  in  forms  similar  to  an  in¬ 
ternal  Kelvin  wave?  To  examine  along  this  line,  we  use 
LZSNFS  outputs  to  eompute  the  daily  average  of  stratifica¬ 
tion  mesh-by-mesh  for  2  June  2007,  and  derive  the  phase 
speed  and  Rossby  radius  of  the  first  baroclinie  mode  around 
Nan  Wan  from  it  (Fig.  5).  Both  quantities  increase  with  the 
water  depth  and  therefore  vary  markedly  in  the  area.  Never¬ 
theless,  relative  to  the  propagation  speed  along  the  axial 


Fig.  4.  Same  as  Fig.  3  except  with  Ihe  removal  of  Coriolis  force  after  2  June  2007  at  0600  GMT. 
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Kig.  5.  7’hc  phase  speed  (a)  and  Rossby  radius  (b)  for  the  first  baroclinic  mode  around  Nan  Wan  on  2  June  2007. 


route  of  cold  anomaly  propagation  as  seen  in  Fig.  3,  Fig.  5 
sugge.sts  a  comparable  phase  speed  of  about  0.5  m  s'*  and  a 
comparable  Rossby  radius  ranging  from  5  to  10  km.  It  is 
difllleult  to  quantify  further  beyond  this  point  because,  in 
addition  to  complexities  mentioned  above,  the  oscillating 
tidal  current  also  modulates  the  propagation. 

For  lack  of  a  better  tema,  we  use  the  term  “internal  Kel¬ 
vin  wave"  loosely  herein  to  include  also  its  nonlinear  and 
variable  bottom  extensions.  In  theory,  the  eft'ect  of  ambient 
current  on  the  internal  Kelvin  wave  propagation  could  be 
nonlinear.  Even  without  the  ambient  current,  nonlinear 
steepening  of  internal  waves  eould  enhance  their  propaga¬ 
tion  speed  up  to  about  1 0  percent.  1  lowever,  the  complex  en¬ 
vironment  around  Nan  Wan.  as  we  illustrated  in  Fig.  4,  pre¬ 
vents  us  from  isolating  the  nonlinear  effect. 

4.  ANOMALOUS  UPWELLINC  IN  JULY  2008 

Figure  6a  shows  the  LZSNFS-produced  sea  level  varia¬ 
tion  (in  blue)  at  station  T  from  27  June  to  8  July  2008.  The 
observed  counterpart  (in  red),  discontinued  for  serv  iee  after 
July  2,  is  also  .superimposed  to  illustrate  the  agreement.  En¬ 
tering  July,  the  tides  begin  to  enter  the  spring  tidal  period.  On 
July  3.  for  example,  the  tidal  range  is  about  1.8  m.  com¬ 
parable  to  the  tidal  range  that  produces  normal  upwclling 
during  spring  tides  (Fig.  2),  but  the  Nan  Wan  upwelling,  as 
we  demonstrate  below,  is  much  stronger. 

Figure  6b  shows  the  observed  bottom  temperature  at 
station  T  for  the  same  period  (in  red)  and  corresponding 
bottom  temperature  produeed  by  LZSNFS  (in  blue).  The 
Novveast/Foreeast  System  reproduces  the  timing  of  each 
temperature  drop  rather  well.  1  lowev  cr.  the  LZSNFS  gener¬ 
ally  underestimates  the  amplitude  of  each  temperature  drop. 
In  particular,  after  each  drop,  the  temperature  often  does  not 


rebound  to  the  observed  level  in  LZSNFS.  Two  reasons  can 
account  for  the  discrepancy.  First,  station  T  is  less  than  I  km 
southwest  of  the  themial  di.sehaige  from  a  nuclear  power 
plant.  Conceivably,  the  thermal  discharge  is  able  to  raise  the 
temperature,  but  its  scale  is  too  small  to  be  included  in 
LZSNFS.  Second,  the  horizontal  resolution  of  LZSNFS 
(2.3  km),  though  fine  enough  from  the  modeling  point  of 
view,  will  likely  generate  horizontal  temperature  diffusion 
that  is  higher  than  reality,  thus  reducing  the  amplitude  of 
temperature  drops  and  rebounds. 

Figure  7  shows  bihourly  currents  and  temperature  at 
25  m  depth  for  2  July  2008.  During  the  time  span,  there  is  a 
strong  subtidal  outllow  from  the  northern  South  China  Sea 
to  the  Philippine  Sea.  When  superimposed  on  the  tidal  cur¬ 
rent,  the  ebb  current  often  exceeds  3  m  s  .  The  outflovv,  tidal 
plus  non-tidal,  is  strong  enough  to  block  the  internal  Kelvin 
wave  propagation,  the  natural  conduit  for  the  dispersal  of 
cold  anomalies.  This  must  be  the  reason  for  the  larger-than- 
normal  temperature  drops. 

At  hour  2,  the  southward  extension  ridge  off  O-Luan-Bi 
lifts  the  eastward  How  to  produce  upwelling.  The  cold  ano¬ 
maly  from  0-Luan-Bi  intensifies,  expands  and  moves  to¬ 
ward  the  central  basin  of  Nan  Wan  from  hour  2  to  hour  6,  and 
dissipates  in  Nan  Wan  at  hour  8.  The  adverse  prevailing  out- 
llovv  is  apparently  too  strong  to  allow  this  cold  anomaly 
moving  to  the  west  coast  of  Taiwan.  At  hour  10,  the  gener¬ 
ally  westward  current  from  east  of  0-Luan-Bi  generates  an¬ 
other  sizeable  cold  anomaly  on  the  lee  (west)  side  of  O- 
Luan-Bi.  It  dissipates  in  central  Nan  Wan  from  hour  10  to 
hour  1 8,  unable  to  move  to  the  west  side  of  Taiwan.  The  cold 
anomaly  generated  off  Mou-Bi-Tou  is  even  more  vulnera¬ 
ble.  At  hour  4,  a  small  region  of  cold  anomaly  first  appears  to 
the  southeast  of  Mou-Bi-Tou.  It  intensifies  and  expands  a 
little  from  hour  4  to  hour  8,  and  dissipates  thereafter.  In 
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general,  cold  anomalies  otTboth  capes  fail  to  disperse  fol¬ 
lowing  the  propagation  of  internal  Kelvin  waves  during  the 
time  span;  blocking  by  the  strong  subtidal  outflow  appears 
to  be  responsible.  Because  of  the  blocking,  temperature 
drops  are  generally  larger  than  nonnal.  The  anomalous  up- 
welling  pattern  as  seen  from  LZSNFS  persists  from  1  July  to 
about  5  July  2008.  This  is  the  period  of  alleged  hypothermia 
in  the  fish  population. 

Figure  8,  similarto  Fig.  5,  shows  the  LZSNFS-produced 
phase  speed  and  Rossby  radius  of  the  ilrst  baroelinic  mode 
around  Nan  Wan  based  on  the  daily  averaged  stratification 
of  2  July  2008.  Relative  to  the  nonnal  upwelling  condition 
on  2  June  2007,  a  slightly  stronger  stratification  on  2  July 
2008  increases  the  phase  speed  and  Rossby  radius  margin¬ 
ally  near  Nan  Wan.  With  such  a  modest  phase  speed  of  about 
0.5  m  s  ’  just  outside  the  perimeter  of  Nan  Wan,  the  strong, 
eastward  subtidal  flow,  as  we  illustrated  in  Fig.  7,  could 
easily  block  the  right-bounded  propagation  of  cold  ano¬ 
malies  in  forms  similar  to  internal  Kelvin  waves. 


5.  SUGGESTED  ORIGIN  OF  ANOMALOUS 
SUBTIDAI.  EASTWARD  FLOW  OFF  NAN  WAN 

Typhoon  Fengshen  moved  through  the  northern  South 
China  Sea  from  22  to  25  June  2008,  leading  us  to  suspect  it 
as  the  main  cause  of  the  anomalous  subtidal  eastward  How 
off  Nan  Wan.  We  now  turn  to  the  larger-scale  EASNFS, 
which  excludes  tides,  to  examine  the  effect  of  this  typhoon. 
Figure  9  contrasts  subtidal  surface  circulation  patterns  be¬ 
tween  normal  upwelling  period  (2  June  2007)  and  anoma¬ 
lous  upwelling  period  (2  July  2008).  The  difference  is  more 
than  subtle.  For  2  June  2007  (Fig.  9a),  the  Kuroshio  more  or 
le.ss  stays  away  from  the  South  China  Sea.  This  is  the  typical 
summerlike  pattern,  in  which  the  Kuroshio  Current  stays 
some di.stance  away  from  Nan  Wan.  For  2  July  2008  (Fig.  9c), 
shortly  after  the  passage  of  typhoon  F'engshen,  the  Kuroshio 
Current  becomes  rather  intrusive;  the  encroachment  pro¬ 
duces  two  anticyclonic  meanders  moving  into  the  northern 
South  China  Sea.  The  northern  meander  produces  strong 


Fig.  6  (a)  The  observed  sea  level  variation  (in  red)  from  27  June  lo  2  July  2(X)8  ai  slaiion  T  and  LZS NFS-produced  counterpart  (in  blue)  from  27  June 
lo  8  July  2(X)8.  (b)  Corresponding  time  series  of b<^)iiom  lemperalure  from  observ  aiion  (red)  and  LZSNFS  (blue).  The  lime  is  in  GMT,  which  is  8  hours 
ahead  of  local  lime. 
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Fig.  7.  Same  as  Fig.  3  except  for  2  July  2008. 
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Kig.  9.  Subtidal  surface  circulation  and  sea  surface  height  (in  color)  for  (a)  normal  upwelling  at  1200  GMT  on  2  June  2007  and  (c)  anomalous 
upwelling  at  l20O(iMT  on  2  July  2008.  prtxiucedby  the  Hast  Asian  Seas  Nowcast/Forecast  System  (KASNFS).  Panels  (b)  and  (d)  show  optimal  inter¬ 
polations  of  corresponding  sea  surface  height  fields  from  all  available,  contemporaneous  altimeter  data. 


.subtkinl  outflow  from  the  northern  South  China  Sea  to  the 
Philippine  Sea  near  Nan  Wan. 

To  the  right  of  Figs.  9a  and  e  are  eontemporaneous 
optimal  interpolations  of  sea  surfaee  height  fields  derived 
from  all  available  satellite  altimeter  measurements,  inelud¬ 
ing  JASON,  GFO,  and  ENVISAT  (Jaeobs  et  al.  2002).  In 
partieular,  this  re.sult  enters  FASNFS  through  data  assimila¬ 
tion,  whieh  slowly  injeets  the  optimal  interpolation  into  the 
dynamieally  active  system  (Chapman  et  al.  2004).  The  simi¬ 
larity  between  EASNFS-derived  and  altimeter-derived  pro¬ 
duets  emerges  as  expeeted  beeause  of  the  data  a.ssimilation. 
The  diserepaney  is  also  understandable  beeause  of  the 
difference  in  methodology.  Similarities  and  dissimilarities 
aside,  one  also  .sees  anomalous  Kuroshio  encroachment  in 
July  2008  from  the  satellite-derived  product  (Fig.  9d). 

The  encroaching  Kuroshio  after  Typhoon  Fengshen  is 
counterintuitive  because  it  is  a  distinctively  winterlike  be¬ 
havior.  Though  different  in  scale,  the  encroachment  is  a 
miniature  form  of  winter  intmsion.  In  winter,  the  prevailing 
northeast  monsoon  moves  waters  from  the  northern  part  to 
the  southern  part  of  the  South  China  Sea.  The  consequent 


low  sea  level  anomaly  in  the  northern  South  China  Sea,  in 
turn,  invites  Kuroshio  to  intrude  (e.g.,  Chem  and  Wang 
2003).  In  other  words,  the  low  sea  level  in  the  northern 
South  China  Sea  is  the  required  preconditioning  for  the 
Kuroshio  intrusion  in  winter.  Extrapolating  the  winter  sce¬ 
nario  further,  did  Typhoon  Fengshen  create  a  similar  but 
weaker,  winter-like  low  sea  level  anomaly  to  invite  the 
Kuroshio  encroachment? 

Figure  10,  again  produced  by  EASNFS,  shows  selected 
snapshots  of  wind  stress  field,  ocean  surface  currents  and 
sea  surface  height  from  1 800  GMT  22  June  to  1 500  GMT  1 
July  2008.  Typhoon  Fengshen  moved  rather  slowly  from 
Luzon  to  southern  China  in  three  days  (from  June  22  to  25). 
A  stationary  typhoon  over  the  northeni  South  China  Sea 
would  propel  waters  away  from  its  center  through  surface 
Ekman  transport  and  generate  a  sea  level  depression  under¬ 
neath  it.  Most  of  the  excess  waters  should  exit  the  Luzon 
Strait,  the  only  major  deep  opening  of  the  South  China  Sea. 
Since  Fengshen  was  moving,  one  expects  a  delayed  re¬ 
sponse  in  its  wake.  A  moving  typhoon  will  also  generate  in¬ 
ertial  currents  radiating  away  from  its  wake  with  a  period  of 
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Fig.  10.  Time  (GMT)  evolution  of  wind  stress  fields  (while  arrows),  ocean  surface  currenls  (black  arrows),  and  ihe  sea  surface  height  dunng  and  after 
ihc  pa.ssagc  of  1  yphoon  l-'cngshcn.  Calibration  arrows  arc  in  Ihc  lop  left  panel 


about  35  hours  in  the  northern  South  China  Sea. 

Figure  10  clearly  .shows  radiating  inertial  currents  from 
the  wake  of  Fengshen.  Of  more  relevance  to  Nan  Wan  is  the 
strong  eastward  subtidal  flow  in  its  vicinity  on  June  26  (mid¬ 
dle  right  panel).  The  eastward  flow  does  not  induce  strong 
iipvvelling  in  Nan  Wan  or  hypothemiia  in  the  fish  population, 
because  June  26  is  in  the  middle  of  neap  tide  period.  After 
Fengshen  landing  on  southern  China,  the  sea  level  depres¬ 
sion  in  the  northern  South  China  Sea  becomes  rather  exten¬ 
sive,  apparently  squeezing  a  substantial  amount  of  surface 
water  away  from  the  South  China  Sea  into  the  Philippine 
Sea.  The  consequent  pressure  deficit  in  the  northern  South 
China  Sea  invites  Kuroshio  encroachment,  inducing  strong, 
eastward  subtidal  flow  near  Nan  Wan  on  July  1  (bottom  right 
panel).  Simply  put,  the  radiating  inertial  current  from  the 
wake  of  Fengshen  first  initiated  the  eastward  subtidal  flow 


near  Nan  Wan,  and  the  Kuroshio  encroachment  continued  it 
later  on 

Obviously,  the  low  sea  level  anomaly  in  the  northern 
South  China  Sea  must  be  close  enough  to  the  Kuroshio  to 
induce  the  Kuroshio  encroachment.  Figure  1 1  (left  panel), 
again  from  EASNFS,  shows  the  zonal  distribution  of  three- 
hourly  sea  surface  height,  averaged  from  20'^N  to  22^N, 
before,  during  and  after  Fengshen  passage.  Fengshen  moved 
through  the  northern  South  China  Sea  from  June  23  to  25. 
Shortly  after,  the  delayed  formation  of  low  sea  level  ano¬ 
maly  expands,  with  a  minor  pause  in  between,  to  the 
Luzon  Strait  ( 1 2 1  °E)  on  July  1 .  Thereafter  the  Kuroshio, 
characterized  by  high  sea  level,  begins  to  encroach  as  an¬ 
ticipated. 

Figure  1 1,  right  panel,  shows  a  concurrent  time  series 
of  sea  surface  height  anomaly  at  21.5'^N  and  120®E  from 
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rig.  1 1 .  Zonal  di.sihbulion  of  sea  surface  height  (averaged  from  20®N  lo  22®N)  before,  during,  and  after  Fengshen  passage,  produced  by  3-hourly  out¬ 
puts  of  EASNFS  (left),  and  saiellite-denved  sea  surface  height  at  2 1 .5°N  and  120®E  from  AVISO  (nght). 


AVISO  NRT  (near  real  time)  product.  AVISO  derives  the 
sea  surface  height  from  near  real  time  altimeter  data  with 
algorithms  that  are  completely  independent  from  the  treat¬ 
ment  of  altimeter  data  in  EASNFS.  It  therefore  provides  a 
somewhat  independent  verification  of  EASNFS's  perfor¬ 
mance.  Specifically,  AVISO  NRT  product  blends  along- 
track  measurements  of  JASON,  ERS-2,  Topex-Poseidon, 
GEO,  and  E VI SAT  with  an  optimal  interpolation.  At  any 
given  time,  the  optimal  interpolation  also  weighs  measure¬ 
ments  in  a  time  w  indow  that  spans  from  43  days  before  to  3 
days  after.  Judging  from  the  result,  AVISO  seems  to  over¬ 
weigh  the  past  and  under-weigh  the  future.  In  consequence, 
the  Kuroshio  encroachment  in  AVISO  occurs  approximately 
one  day  earlier  than  in  EASNFS.  The  slight  discrepancy, 
however,  does  not  upset  the  dynamic  process  we  intend  to 
establish  herein. 

fo  quantify  further.  Fig.  12  .shows  the  3-hoiirly  volume 
transport  across  the  entire  cross-section  of  Luzon  Strait  from 
15  June  to  10  July  2008,  produced  by  the  EASNFS.  Also 


shown  is  a  smoother  curve  after  applying  a  2-day  low-pa.ss 
filter.  In  computing  them,  the  depth  average  is  not  ideal  be¬ 
cause  it  also  contains  contributions  from  greater  depths  that 
are  not  relevant  to  the  Kuroshio  encroachment.  It  neverthe¬ 
less  tells  a  little  more.  Leaving  high-frequency  fluctuations 
aside,  the  inflow  from  the  Philippine  Sea  to  the  South  China 
Sea  generally  increases  prior  to  the  arrival  of  Typhoon 
Fengshen  (June  22).  After  Fengshen  impacts  fully  on  the 
northern  South  China  Sea  (June  24),  the  inflow  from  the 
Philippine  Sea  decreases  continuously  and  eventually  be¬ 
comes  an  outflow  around  June  26.  Apparently,  Fengshen 
creates  a  radiating  inertial  current  in  its  wake  and  propels 
northern  South  China  waters  out  of  the  Luzon  Strait  in  this 
time  span.  The  outflow  reverses  to  become  a  pronounced 
inflow  on  June  29;  this  preconditions  the  Kuroshio  encro¬ 
achment  in  response  to  the  low'  sea  level  of  the  northern 
South  China  Sea.  The  inflow  reverses  back  to  a  weak  out¬ 
flow  shortly  before  July  1  and  gradually  becomes  a  modest 
inflow  again  thereafter. 
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Fig.  12.  Three-hourly  volume  transpi^rt  across  the  Luzon  Strait  ( 1 20.8® K)  before,  dunng  and  alter  the  passiige  oftyphcKin  Fengshen,  produced  by  F^isl 
Asian  Seas  Nowcasty'Forcca.si  System  (FASNKS).  Also  shown  is  the  low-frequency  response  after  applying  a  2-day  low-pass  filter.  Positive  (nega¬ 
tive)  transports  are  eastward  (westward). 


Figure  1 3  shows  power  spectrum  densities  for  the  Litzon 
Strait  transport.  Tlic  red  speetrum,  covering  the  i^eriod  from 
2008/06/01  00  GMT  to  2008/07/18  18  GMT,  contains  a 
substantial,  relatively  quiescent,  pre-typhoon  time  span. 
Two  peaks  stand  out  eounterintuitively  near  diurnal  and  se¬ 
midiurnal  periods,  even  though  the  EASNFS  does  not  simu¬ 
late  tides.  The  near-diurnal  oscillations  arise  from  resonance 
(quarter-wave  seiche)  that  involves  the  entire  South  China 
Sea  (Jan  et  al.  2007).  In  the  Luzon  Strait,  the  peak-to-peak 
zonal  distance  between  the  two  north-south  running  ridges 
(Fig.  la)  is  about  120  km.  It  takes  about  12  hours  for  long 
(non-dispersive)  internal  waves  to  travel  from  one  ridge  to 
the  other  (Chao  ct  al.  2007),  producing  the  near-semidiurnal 
.spectral  peak. 

Restricting  attention  to  the  post-typhoon  period,  the 
blue  spectrum  in  Fig.  13  covers  the  time  span  from  0000 
GM  r  2008/06/25  to  1 800  GMT  2008/07/09.  In  addition  to 
the  two  peaks  near  diunial  and  .semidiurnal  periods,  a  third 
peak  near  the  inertial  period  (35  hours)  emerges.  The  iner¬ 
tial  current  radiating  aw'ay  from  the  w  ake  of  the  typhoon- 
covered  area  produces  these  oseillations  in  the  Luzon  Strait. 


6.  CONCLUSIONS 

Nan  Wan  upwelling  intensifies  during  the  spring  tide. 
Both  capes  and  their  southward  extension  ridges  are  re¬ 
sponsible  for  the  upwelling.  Using  NRL  Luzon  Strait  Ocean 
Now'cast^Forccast  System,  we  have  identified  and  con¬ 
trasted  two  modes  of  dispersal  for  the  upw  elled  cold  ano¬ 
malies.  Under  normal  conditions,  cold-water  masses  pro¬ 
pagate  away  from  Nan  Wan  in  forms  similar  to  internal 
Kelvin  waves.  Since  this  is  the  natural  w  ay  of  dispersion. 


our  results,  although  limited  so  far,  indicate  that  the  eold 
anomalies  arc  less  likely  to  ovcr-accumulate  in  Nan  Wan 
under  this  normal  mode  of  dispersion.  For  early  July  of 
2008,  the  Now'cast/Forecast  System  suggested  the  existence 
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Fig.  13.  Power  spectrum  densities  for  Ihe  3-hourly  Luzon  Strait  trans¬ 
port.  The  red  line  u.ses  a  time  series  from  0000  GMT  2008/06/01  to 
1 800  GMT  2008/07/ 1 8,  including  a  substantial  quiesceni  period  before 
the  arrival  of  Typhoon  Fengshen.  The  blue  line  uses  a  shorter,  posl- 
lyphoon  time  scries  from  OOOO  GMT  2008/06/25  to  1800  (iMT 
2008/07/09,  accentuating  the  |x>st-typhoon  inertial  oscillations. 
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of  an  anomalous,  subtidal  and  essentially  easm^ard  flow  near 
Nan  Wan,  which,  when  superimposed  on  the  ebb  tide,  blocks 
the  cold  water  dispersion  in  the  propagation  direction  of  in¬ 
ternal  Kelvin  waves.  Wc  suggested  Typhoon  Fengshen  as 
the  main  cause  of  the  anomalous  subtidal  flow.  Under  this 
anomalous  condition,  the  Nowcast/Forecast  System  sug¬ 
gested  over-accumulation  of  cold  anomalies  in  Nan  Wan. 
Not  surprisingly,  the  alleged  hypothcmiia  in  the  fish  popula¬ 
tion  occurred  during  this  period. 

Our  results,  though  a  step  fonvard  in  understanding  the 
Nan  Wan  upwelling,  also  open  a  series  of  questions  waiting 
for  answers.  For  example,  are  there  still  other  modes  of  up¬ 
welling  in  Nan  Wan?  Will  typhoons  with  substantially  dif¬ 
ferent  tracks  also  induce  the  anomalous  upwelling  in  Nan 
Wan?  During  the  least  observed  but  extensive  Kuroshio  in- 
tnision  period  in  winter,  is  Nan  Wan  upwelling  mostly  as 
anomalous  as  in  July  2008?  Our  investigations  along  this 
line  are  ongoing. 
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